Energy scales and magnetoresistance at a quantum critical point 
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The magnetoresistance (MR) of CeCoIns is notably different from that in many conventional met- 
als. We show that a pronounced crossover from negative to positive MR at elevated temperatures 
and fixed magnetic fields is determined by the scaling behavior of quasiparticle effective mass. At a 
quantum critical point (QCP) this dependence generates kinks (crossover points from fast to slow 
growth) in thermodynamic characteristics (like specific heat, magnetization etc) at some tempera- 
tures when a strongly correlated electron system transits from the magnetic field induced Landau 
Fermi liquid (LFL) regime to the non-Fermi liquid (NFL) one taking place at rising temperatures. 
We show that the above kink-like peculiarity separates two distinct energy scales in QCP vicinity - 
low temperature LFL scale and high temperature one related to NFL regime. Our comprehensive 
theoretical analysis of experimental data permits to reveal for the first time new MR and kinks 
scaling behavior as well as to identify the physical reasons for above energy scales. 
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INTRODUCTION 

An explanation of rich and striking behavior of 
strongly correlated electron system in heavy fermion 
(HF) metals is, as years before, among the main problems 
of condensed matter physics. One of the most interesting 
and puzzling issues in the research of HF metals is their 
anomalous normal-state transport properties. Measure- 
ments of magnetoresistance (MR) on CeCoIns P, Q have 
shown that it is notably different from ordinary weak- 
field orbital MR described by Kohler's rule which holds 
in many conventional metals, see e.g. Q. At fixed mag- 
netic fields B, MR of CeCoIns exhibits a crossover from 
negative (low temperatures) to positive (hight tempera- 
tures) one at temperature growth [l|, Q ■ This crossover is 
hard to explain within both conventional Fermi liquid ap- 
proach for metals and in terms of Kondo systems [4] . To 
explain this effect, it has been assumed that the crossover 
can be attributed to some distinct energy scales revealed 
by kinks (crossover points from fast to slow growth) in 
thermodynamic characteristics (like specific heat, magne- 
tization etc) and leading to a change of spin fluctuations 
character with increasing of the applied magnetic field 
strength 1^, \4, \^ S, i, i] . 

Here we investigate the NFL-LFL transition region 
(we call it below crossover region), where MR changes 
its sign. The modified Kohler's rule (MR versus tan- 
gent of Hall angle) have been utilized to describe MR 
data d, In this region, both Kohler's rule and its 
modified version do not work. In Landau Fermi liquid 
(LFL) regime, the quasiparticles were observed in mea- 
surements of transport properties in CeCoIns [loj . An 
analysis of above thermodynamic quantities shows that 
quasiparticles exist in both LFL and crossover regimes 



when strongly correlated Fermi swtems like HF met- 
als or two-dimensional (2D) ^He [nl, [H, [H, Q [H, [3 
transit from LFL to NFL behavior. It is of crucial im- 
portance to verify whether quasiparticles with effective 
mass Af * still exist and determine the transport prop- 
erties and energy scales in HF metals in crossover re- 
gion. On the other hand, even early measurements on 
HF metals gave evidences in favor of the quasiparticles 
existence. For example, the application of magnetic field 
B restores LFL behavior of HF metals which demon- 
strate NFL properties in the absence of the field . In 
that case the empirical Kadowaki- Woods (KW) ratio K 
is conserved, K = A{B)/j^{B) cx A{B)/x'^{B) = const 
C/T, 
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where 70 = C/T, C is a heat capacity, x is 
a magnetic susceptibility and A{B) is a coefficient de- 
termining the temperature dependence of the resistivity 
p = po + A{B)T'^ . Here po is the residual resistance. 
The observed conservation of K can be hardly interpreted 
within scenarios when quasiparticles are suppressed, for 
there is no reason to expect that 70 (S), x{T), A{B) and 
other transport and thermodynamic quantities like ther- 
mal expansion coefficient a{B) are affected by the fluc- 
tuations or localization in a correlated fashion. As we 
will see below, the MR measurements in the crossover 
region can present indicative data on the quasiparticles 
availability. Such MR measurements were carried out 
in CeCoIns when the system transits from LFL to NFL 
regime at elevated temperatures and fixed magnetic fields 

m 

In this Letter, we analyze MR of CeCoIns and show 
that the crossover from negative to positive MR at ele- 
vated temperatures and fixed magnetic fields can be well 
captured utilizing fermion condensation quantum phase 
transition (FCQPT) concept based on the quasiparti- 
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cles paradigm 1^, |20|, |2l|, [2^. We demonstrate that 
the crossover is regulated by the universal behavior of 
the effective mass M*{B,T) observed in many HF met- 
als. It is exhibited by M*{B,T) when HF metal transits 
from LFL regime (induced by the application of magnetic 
field) to NFL one taking place at rising temperatures. 
The above behavior of the effective mass also generates 
kinks (crossover points from fast to slow growth at el- 
evated temperatures) in thermodynamic characteristics 
(like specific heat, magnetization etc). We show that the 
above kink-like peculiarity separates two distinct energy 
scales - low temperature LFL scale and high tempera- 
ture one related to NFL regime. Our calculations of MR 
are in good agreement with observations and allow us to 
reveal new scaling behavior of both MR and the kinks. 



SCALING BEHAVIOR OF THE KINKS 

To study universal low temperature features of HF 
metals, we use the model of homogeneous heavy-fermion 
liquid with the effective mass M*{T,B,x), where x = 
Pf/Stt^ is a number density and pf is Fermi momentum 
[23|. This model permits to avoid complications associ- 
ated with the crystalline anisotropy of solids [l^. We 
first outline the case when at T ^ the heavy-electron 
liquid behaves as LFL and is brought to the LFL side 
of FCQPT by tuning of a control parameter like x. At 
elevated temperatures the system transits to the NFL 
state. The dependence Af*(T, x) is governed by Landau 
equation j23j 
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where n(p, T, x) is the distribution function of quasipar- 
ticles and -F(pi?,p) is Landau interaction amplitude, M 
is a free electron mass. At T = 0, eq. (P) reads [2^ 
M*/M = 1/(1 - NaF^{pF,PF)/S)- Here A^o is the den- 
sity of states of a free electron gas, F^{pf,Pf) is the 
p-wave component of Landau interaction amplitude F. 
Taking into account that x — pI^/Stt^, we rewrite the am- 
plitude as F^{pf,Pf) = F^{x). When at some critical 
point X = xpc, F^{x) achieves certain threshold value, 
the denominator tends to zero and the system under- 
goes FCQPT related to divergency of the effective mass 
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Equation ([2]) is valid in both 3D and 2D cases, while the 
values of factors A and B depend on the dimensionality. 
The approximate solution of Eq. (jlj is of the form [l4| 
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where Ai > and A2 < are constants of order unity, 
(3 = TM*{T)/pI and /(O) - ^^(xfc)- K follows from 
Eq. ^ that the effective mass M* as a function of T 
and X reveals three different regimes at growing temper- 
ature. At the lowest temperatures we have LFL regime 
with M*(T, x) ~ M*{x) + aT^ with a < since Ai > 0. 
The effective mass as a function of T decays down to a 
minimum and after grows, reaching its maximum M^^j 
at some temperature Tm(x) then subsequently dimin- 
ishing as T"^/'^ [III, Moreover, the closer is the 
number density x to its threshold value xfci the higher 
is the rate of the growth. The peak value M^j- grows 
also, but the maximum temperature Tm lowers. Near 
this temperature the last "traces" of LFL regime disap- 
pear, manifesting themselves in the divergence of above 
low-temperature series and substantial growth of M*{x). 
Numerical calculations based on Eqs. ([T]) and ([3]) show 
that at rising temperatures T > T1/2 {T1/2 is a charac- 
teristic temperature determining the validity of regime 
(III), see Ref. [l3| for details) the linear term (x (3 gives 
the main contribution and leads to new regime when Eq. 
© reads M/M*{T) cx (3 yielding 



M*{T) cx T 
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We remark that Eq. (HI ensures that at T > T1/2 the 
resistivity behaves as p{T) cx T [12]. Near the critical 
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FIG. 1: Schematic phase diagram of the systems under consid- 
eration. Control parameter C, represents number density (or 
doping) s, magnetic field B, pressure P etc. i^fc denotes the 
point of effective mass divergence. SC, FM, AFM denote the 
superconducting, ferromagnetic and antiferromagnetic states, 
respectively. Thp vertical arrow shows LFL-NFL transitions 
at rising temperatures and fixed C,. Inset shows a schematic 
plot of the normalized effective mass M'^ = M* [T /Tm) / M^i 
{Mil is its maximal value at T = Tm) versus the normalized 
temperature Tjv = T/Tm- Several regions are shown. First 
goes the LFL regime {M1;{Tn) ~ const) at Tn <S 1, then 
the transition regime (the hatched area) where Mj^ reaches 
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its maximum. At elevated temperatures T 
followed by T"^''^ behavior, see Eq. (gj. 
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regime occurs 
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point Xfc {M/M*{x — > xpc) 0), the behavior of the 
effective mass changes dramatically since the first term 
in the right-hand side of Eq. ^ vanishes so that the 
second term becomes dominant. As a result, we can no 
more measure the mass M* in units of M (as M/M*{x 
Xfc) 0) and we have to measure M* in units of Af|^ 
and T in units of T^/. Latter scales can be viewed as 
natural ones. 

The schematic phase diagram of HF liquid is reported 
Fig. [T] The control parameter C can be pressure P, mag- 
netic field B, or doping (density) x. At C = Cfc, FCQPT 
takes place leading to a strongly degenerated state. This 
state is captured by the superconducting (SC), ferromag- 
netic (FM), antiferromagnetic (AFM) etc. states lifting 
the degeneracy [13] ■ The variation of C drives the sys- 
tem from NFL region to LFL one. For example, in the 
case of magnetic field B, (pc = -Bco, where Bco is a crit- 
ical magnetic field, such that at _B > Bco the system is 
driven towards its LFL regime. Below we consider the 
case with ( > ^pc when the system is on the LFL side 
of FCQPT. The inset demonstrates the behavior of the 
normalized effective mass = M* /M^ versus normal- 
ized temperature Tn = T/Tm- Both T-^/^ and T'^/^ 
regimes are marked as NFL ones since the effective mass 
depends strongly on temperature. The temperature re- 
gion T ~ Tm signifies the crossover between the LFL 
regime with almost constant effective mass and NFL be- 
havior, given by T~'^/^ dependence. Thus temperatures 
T ~ Tm can be regarded as the crossover region between 
LFL and NFL regimes. 

It turns out that M* (T, x) in the entire T <Ti/2 range 
can be well approximated by a simple universal inter- 
polating function [ll|, [l^, ■ The interpolation occurs 
between the LFL (M* cx T^) and NFL (M* cx T'^'^, 
see Eq. (j4|)) regimes thus describing the above crossover. 
Introducing the dimensionless variable y = T/v = TjTu-, 
we obtain the desired expression 
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(5) 



Here M^{y) is the normalized effective mass, ci and C2 
are parameters, obtained from the condition of best fit to 
experiment. To correct the behavior of M^(j/) at rising 
temperatures M* ^ T~^/^, we add a term to Eq. ([5]) 
and obtain 
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where C3 is a parameter. The last term on the right hand 
side of Eq. ([6]) makes satisfy Eq. (|4]) at temperatures 
T/Tm > 2. 

At small magnetic fields B (that means that Zeeman 
splitting is small), the effective mass does not depend on 
spin variable and B enters Eq. ([T]) as Bllb/T {fXB is 
Bohr magneton) making Tm oc B^b [HI, The 



application of magnetic field restores the LFL behavior 
and at T < Tm the effective mass depends on _B as 11 
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M*{B) cx {B-B,o) 
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(7) 



Note that in some cases Bco — 0. For example, the HF 
metal CeRu2Si2is characterized by Bco — and shows 
neither evidence of the magnetic ordering or supercon- 
ductivity nor the LFL behavior down to the lowest tem- 
peratures [24| . In our simple model Bcq is taken as a 
parameter. We conclude that under the application of 
magnetic field the variable 



y = T/Tm oc 



T 



(8) 



f^siB — Bco) 

remains the same and the normalized effective mass is 
again governed by Eqs. ([5]) and ^ which are the final 
result of our analytical calculations. We note that the 
obtained results are in agreement with numerical calcu- 
lations [lll,[ll. 

The normalized effective mass M^{y) can be ex- 
tracted from experiments on HF metals. For example, 
M*{T,B) cx C(T)/T oc S{T)/T cx xac[T), where S{T) 
is the entropy, C(T) is the specific heat and xac{T) is 
ac magnetic susceptibility [ll|i [I2| ■ If the corresponding 
measurements are carried out at fixed magnetic field B 
(or at fixed x and B) then, as it seen from Fig. [1] the 
effective mass reaches the maximum at some tempera- 
ture Tm- Upon normalizing both the effective mass by 
its peak value Ml,j at each field B and the temperature 
by Tm, we observe that all the curves merge into a sin- 
gle one, given by Eqs. (O and Q thus demonstrating a 
scaling behavior. 

To verify Eq. ([3]), we use measurements of xac{T) 
in CeRu2Si2 at B = 0.02 mT at which this HF metal 



demonstrates the NFL behavior [2^ . It is seen from Fig. 
[5] that Eq. ^ gives good description of the facts in the 
extremely wide range of temperatures. The inset to Fig. 
[5] exhibits a fit for M^{y) extracted from measurements 
oiXAciT) at different magnetic fields, clearly indicating 
that the function given by Eq. ([5]) represents a good 
approximation for M^{y) when the system transits from 
the LFL regime to NFL one. 

M^{y) extracted from the entropy S[T)/T and mag- 
netization M measurements on the "^He film [5^ at dif- 
ferent densities x is reported in the left panel of Fig. [S] 
In the same panel, the data extracted from the heat ca- 
pacity of the ferromagnet CePdo.2R-ho.8 [1^ and the AC 
magnetic susceptibility of the paramagnet CeRu2Si2 24 1 
are plotted for different magnetic fields. It is seen that 
the universal behavior of the normalized effective mass 
given by Eq. ([5]) and shown by the solid curve is in ac- 
cord with the experimental facts. All 2D ■^He substances 
are located at C > Cfc (see Fig. [TJ, where the system 
progressively disrupts its LFL behavior at elevated tem- 
peratures. In that case the control parameter, driving 
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FIG. 2: Temperature dependence of the ac susceptibility xac 
for CeRu2Si2. The solid curve is a fit for the data shown 
by the triangles at _B = 0.02 mT and represented by the 
function x{T) = a/y/T given by Eq. Q with a being a fitting 
parameter. Inset shows the normalized effective mass versus 
normalized temperature y extracted from xac measured at 
difi^erent fields as indicated in the inset [23 ]. The solid curve 
traces the universal behavior of M%{y) determined by Eq. 
((5]). Parameters ci and C2 are adjusted to fit the average 
behavior of the normalized effective mass M^{y). 



the system towards its quantum critical point (QCP) is 
merely the number density x. It is seen that the behav- 
ior of M^{y), extracted from S{T)/T and magnetization 
M of 2D '^He looks very much like that of 3D HF com- 
pounds. In the right panel of Fig. [31 the normalized data 
on C{y), S{y), yxiu) and M{y) + yxiv) extracted from 
data collected on CePdi^xRh^ [26^ , ^He [^H, CeRusSia 
[2^ . CeCoIns 27] and YbRu2Si2 [Tj] respectively are pre- 
sented. Note that in the case of YbRu2Si2, the variable 
y = {B — Bco)^iB /Tm can be viewed as effective normal- 
ized temperature. As seen from Eq. (O, this representa- 
tion of the variable y is correct when the temperature is 
a fixed parameter. 

It is seen from the right panel of Fig. [3] that all the 
data exhibit the kink (shown by arrow) at y > 1 taking 
place as soon as the system enters the transition region 
from the LFL regime to the NFL one and corresponding 
to the temperatures where the vertical arrow in Fig. [T] 
crosses the solid line. It is also seen that the low tem- 
perature LFL scale of the thermodynamic functions (as 
a function of y) is characterized by the fast growth and 
the high temperature one related to the NFL behavior 
is characterized by the slow growth. As a result, we can 
identify the energy scales near QCP, discovered in Ref. 
[3]: the thermodynamic characteristics exhibit the kinks 
(crossover points from the fast to slow growth at elevated 
temperatures) which separate the low temperature LFL 
scale and high temperature one related to NFL regime. 



FIG. 3: The left panel. The normalized effective mass M% 
versus the normalized temperature y = T/Tm- The depen- 
dence M^ijj) is extracted from measurements of S{T)/T and 
magnetization M on 2D ^He [25ll) . from ac susceptibility 
XAc{T) collected on CeRu2Si2 [3 and from C{T)/T col- 
lected on CePdi_xRhx |2^. The data are collected for differ- 
ent densities and magnetic fields shown in the left bottom cor- 
ner. The solid curve traces the universal behavior of the nor- 
malized effective mass determined by Eq. ©. Parameters ci 
and C2 are adjusted for xn{Tn, B) at B = 0.94 mT. The right 
panel. The normalized specific heat C{y) of CePdi_xRhx at 
different magnetic fields _B, normalized entropy S{y) of ^He 
at different number densities x, and the normalized yx(y) at 
B — 0.94 mT versus normalized temperature y are shown. 
The upright triangles depict the normalized 'average' magne- 
tization M + Bx collected on YbRu2Si2 0- The kink (shown 
by the arrow) in all the data is clearly seen in the transi- 
tion region y > 1. The solid curve represents yM^ijj) with 
parameters ci and C2 adjusted for magnetic susceptibility of 
CeRu2Si2 at B = 0.94 mT. 



SCALING BEHAVIOR OF THE 
MAGNETORESISTANCE 

By definition, MR is given by 

p(S,T)-p(0,T) 



Pmr{B, T) 



(9) 



We apply Eq. ^ to study MR of strongly correlated elec- 
tron liquid versus temperature T as a function of mag- 
netic field B. The resistivity p{B,T) is 

piB,T)=po + ApiB,T) + ApLiB,T), (10) 

where po is a residual resistance, Ap — ciAT^, ci is a 
constant, A is a coefficient determining the temperature 
dependence of the resistivity p = po + AT'^ . The classical 
contribution Apl{B,T) to MR due to orbital motion of 
carriers induced by the Lorentz force obeys the Kohler's 
rule We note that Apl{B) <^ p(0, T) as it is assumed 
in the weak-field approximation. To calculate A, we use 
the quantities 70 = C/T oc M* and/or x oc M* as well as 
employ the fact that Kadowaki- Woods ratio K = A/jq oc 
A/x'^ = const 113]. As a result, we obtain A cx (M*)^ 
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T3,[ii[li, so that Ap{B,T) = c{M* {B ,T) fT'^ and c is 
a constant. Suppose that the temperature is not very low, 
so that po < Ap{B = 0,T)^nd B > B^o- Substituting 
(fTUI into dH), we find that ^ 



ApL 



cT 



{M*{B,T)Y - (M*(0,r))2 
P(0,T) 



(11) 



Consider the quahtative behavior of MR described by 
Eq. pip as a function of B at a certain temperature 
T ~ Tq. In weak magnetic fields, when Tq > T1/2 and the 
system exhibits NFL regime (see Fig. [IJ, the main con- 
tribution to MR is made by the term Apl{B), because 
the effective mass is independent of the applied magnetic 
field. Hence, \M*{B,T) ~ M*{0,T)\/M*{0,T) < 1 and 
the leading contribution is made by Apl{B). As a result, 
MR is an increasing function of B. When B becomes so 
high that Tm{B) ~ pb{B — -Bco) ~ To, the difference 
{M*{B,T) - M*(0,r)) becomes negative and MR as a 
function of B reaches its maximal value at Tm{B) ^ Tq 
when the kink occurs, see the right panel of Fig. [31 At 
further increase of magnetic field, when Tm{B) > Tq, the 
effective mass M*{B,T) becomes a decreasing function 
of B, as follows from Eq. As B increases. 



(M*(B,T) -M*(0,r)) 
M*{0,T) 



(12) 



and the magnetoresistance, being a decreasing function 
of B, is negative. 

Now we study the behavior of MR as a function of 
T at fixed value Bq of magnetic field. At low tempera- 
tures T < Tm{Bq), it follows from Eqs. © and ^ that 
M*{Bo,T)/M*{0,T) < 1, and it is seen from Eq. ^ 

that pmriBo,T) 1, because Apl{Bq,T)/ p{Q,T) < 

1. We note that Bq must be relatively high to guar- 
antee that M*{Bq,T)/M*{0,T) < 1. As the tem- 
perature increases, MR increases, remaining negative. 
At T ~ Tm{Bo), MR is approximately zero, because 
p{Bq,T) ~ p{Q,T) at this point. This allows us to con- 
clude that the change of the temperature dependence of 
resistivity p{Bo, T) from quadratic to linear manifests it- 
self in the transition from negative to positive MR. One 
can also say that the transition takes place when the kink 
occurs (as shown by the arrow in the right panel of Fig. 
[3]) and the system goes from the LFL behavior to the 
NFL one. At T > Tm{Bq), the leading contribution to 
MR is made by Apl{Bo,T) and MR reaches its maxi- 
mum. At Tm{Bo) <^ T, MR is a decreasing function of 
the temperature, because 



\M*{B,T) - M*{0,T)\ 
M*(0,T) 



« 1, 



(13) 



and pmr{Bo,T) <^ 1. Both transitions (from positive to 
negative MR with increasing B at fixed temperature T 
and from negative to positive MR with increasing T at 




Normalized temperature 

FIG. 4: The normalized magnetoresistance R%(y) given by 
Eq. (|15|l versus normalized temperature y = T/TRm. R'^{y) 
was extracted from MR shown in Fig. [6] and collected on 
CeCoIns at fixed magnetic fields B [3) listed in the right upper 
corner. The starred line represents our calculations based on 
Eqs. ([5]) and (|15|l with the parameters extracted from ac 
susceptibility of CeRu2Si2 (see the caption to Fig. [2]). The 
solid line displays our calculations based on Eqs. (O and 
(|15p : only one parameter was used to fit the data, while the 
other were extracted from the ac susceptibility measured on 
CeRu2Si2. 



fixed B value) have been detected in measurements of 
the resistivity of CeCoIns in a magnetic field [ij . 

Let us turn to quantitative analysis of MR. As it was 
mentioned above, we can safely assume that the classical 
contribution Apl{B,T) to MR is small as compared to 
Ap{B,T). Omission of Apl{B,T) allows us to make 
our analysis and results transparent and simple since the 
behavior of Ap]^{Bq,T) is not known in the case of HF 
metals. Consider the ratio Rp = p{B,T)/ p{0,T) and 
assume for a while that the residual resistance po is small 
in comparison with the temperature dependent terms. 
Taking into account Eq. (|10|) and p(0, T) oc T, we obtain 
fromEq. (HH) 



R" = Pn 



1 = 



P{B,T) 

p(o,r) 



cx T{M*{B,T)f 



(14) 



It follows from Eqs. ([5]) and that the ratio Rp reaches 
its maximal value i?^^ at some temperature ~ ?a/- 
If the ratio is measured in units of its maximal value i?^^ 
and T is measured in units of T^m ^ Tm then it is seen 
from Eqs. (O, ^ and that the normahzed MR 



RP{B,T) 



(15) 



becomes a universal function of the only variable y = 
T/Tum- To verify Eq. ^5]) . we use MR obtained in mea- 
surements on CeCoIns, see Fig. 1(b) of Ref. [l|. The re- 
sults of the normalization procedure of MR are reported 
in Fig. m It is clearly seen that the data collapse into 
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the same curve, indicating that the normahzed magne- 
toresistance i?^ well obeys the scaling behavior given by 
Eq. This scaling behavior obtained directly from 

the experimental facts is a vivid evidence that MR be- 
havior is predominantly governed by the effective mass 
M*(B,T). 

Now wc are in position to calculate R'^ii/) given by 



Eq. (HH). Using Eq. ^ to parameterize M^{y), we 
extract parameters ci and C2 from measurements of the 
magnetic ac susceptibility x on CeRu2Si2 [24| and apply 
Eq. (|15p to calculate the normalized ratio. It is seen 
that the calculations shown by the starred line in Fig. 
S] start to deviate from experimental points at elevated 
temperatures. To improve the coincidence, we employ 
Eq. ^ which describes the behavior of the effective mass 
at elevated temperatures in accord with Eq. ^ and 
ensures that at these temperatures the resistance behaves 
as p{T) cx T. In Fig. |1 the fit of i?^(y) by Eq. ^ 
is shown by the solid line. Constant C3 is taken as a 
fitting parameter, while the other were extracted from 
ac susceptibility of CeRu2Si2 as described in the caption 
to Fig. H 
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FIG. 5: The peak temperatures TRm (squares) and the peak 
values -Rmax (triangles) versus magnetic field B extracted from 
measurements of MR l] . The solid lines represent our calcu- 
lations based on Eqs. (|16p and (|17|l . 

Before discussing the magnetoresistance pmr{B,T) 
given by Eq. ([9]), we consider the magnetic field depen- 
dencies of both the MR peak value i?max(-B) and corre- 
sponding peak temperature T^iraiB). It is possible to use 
Eq. pi)) which relates the position and value of the peak 
with the function M*{B, T). Since Tptm oc psB, B enters 
Eq. p4)) only as tuning parameter of QCP, as both Ap^ 
and Pa were omitted. At i? ^ i?co and T < T^ni{B), this 
omission is not correct since Ap^ and pa become com- 
parable with Ap{B,T). Therefore, both i?max(S) and 
TnmiB) are not characterized by any critical field, being 
a continuous function at the quantum critical filed Bco, in 
contrast to M*{B,T) which peak value diverges and the 
peak temperature tends to zero at Bco as it follows from 
Eqs. ([7]) and Thus, we have to take into account 
Api(_B, T) and po which prevent T^jniB) from vanishing 
and make i?max(-B) finite at _B ^ Bco- As a result, we 



have to replace Bco by some effective field B^ff < Bco 
and take B^ff as a parameter which imitates the con- 
tributions coming from both Apl{B,T) and po- Upon 
modifying Eq. (|I4p by taking into account Apl{B,T) 
and Po, we obtain 



TRm(B) C 
Rmax{B) 



h{B-B,ff), 
^ h2{B-B,ff ) 



-1/3 



63(i?-i3e//)-l + l 



(16) 

(17) 



Here 61, 62, 63 and -Bg// are fitting parameters. It is 
pertinent to note that while deriving Eq. (|17p we use 



Eq. mD with substitution {B - B^ff) for T. Then, 
Eqs. ([HI) and HT]) are not vahd at _B < Bcq. In Fig. 
[5l we show the field dependence of both Tr,,! and i^max, 
extracted from measurements of MR [1]. It is seen that 
both Tftni and i?max are well described by Eqs. and 
P7)) with Be/ / =3.8 T. We note that this value of B^f / is 
in good agreement with observations obtained from the 
B — T phase diagram of CeCoIus, see the position of the 
MR maximum shown by the filled circles in Fig. 3 of Ref. 



a. 




T[K] 

FIG. 6: MR versus temperature T as a function of magnetic 
field B. The experimental data on MR were collected on 
CeCoIus at fixed magnetic field B [J] shown in the right bot- 
tom corner of the Figure. The solid lines represent our calcu- 
lations, Eq. ([5} is used to fit the effective mass entering Eq. 

To calculate Pmr{B,T), we apply Eq. ^5)1 to describe 
its universal behavior, Eq. (O for the effective mass along 
with Eqs. pB)) and P?|) for MR parameters. Figure [H 
shows the calculated MR versus temperature as a func- 
tion of magnetic field B together with the experimental 
points from Ref. [ij. We recall that the contributions 
coming from Apl{B ,T) and po were omitted. As seen 
from Fig. [51 our description of experiment is pretty good. 

SUMMARY 

Our comprehensive theoretical study of MR shows that 
it is (similar to other thermodynamic characteristics like 
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magnetic susceptibility, specific heat etc) governed by the 
scahng behavior of the quasiparticle effective mass. The 
crossover from negative to positive MR occurs at elevated 
temperatures and fixed magnetic fields when the system 
transits from the LFL behavior to NFL one and can be 
well captured by this scaling behavior. This behavior 
permits to identify the energy scales near QCP, discov- 
ered in Ref. Namely, the thermodynamic charac- 
teristics (like specific heat, magnetization etc) consist of 
the low temperature LFL scale characterized by the fast 
growth and the high temperature one related to the NFL 
behavior and characterized by the slow growth. These 
scales are separated by the kinks in the transition region. 
Obtained theoretical results are in good agreement with 
experimental facts and allow us to reveal for the first 
time a new scaling behavior of both magnetoresistance 
and kinks separating the different energy scales. 
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